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^ ■ We consider the problem of the definition of an effective temperature via 

the long-time limit of the fluctuation-dissipation ratio X°° after a quench 
from the disordered state to the critical point of an O(N) model with dissi- 
— L ■ pative dynamics. The scaling forms of the response and correlation functions 

Ph ■ of a generic observable 0(t) are derived from the solutions of the correspond- 

ir . ing Renormalization Group equations. We show that within the Gaussian 

approximation all the local observables have the same Xq , allowing for a 
£\j \ definition of a unique effective temperature. This is no longer the case when 

fluctuations are taken into account beyond that approximation, as shown by 
qq ■ a computation up to the first order in the e-expansion for two quadratic ob- 

servables. This implies that, contrarily to what often conjectured, a unique 
effective temperature can not be defined for this class of models. 



O 

o 

03 



-a 
c 
o 
o 



X 
S3 



PACS Numbers: 05.10.Cc, 05.10.Gg, 64.60.Ht, 75.40.Gb 



Typeset using REVTeX 



I. INTRODUCTION 

The non-equilibrium dynamics of physical systems is one of the most challenging prob- 
lems in statistical mechanics. Equilibrium statistical mechanics has been probably one of the 
most important achievements during the last century. On the other hand, in nature equi- 
librium is more an exception rather than a rule. In view of that many efforts are currently 
aiming at achieving a coherent theoretical picture of non-equilibrium phenomena: Indeed 
many real systems persist out-of-equilibrium practically forever. One example is naturally 
provided by slow-relaxing systems, such as structural glasses and spin-glasses, whose equi- 
libration times sometimes exceed any reasonable experimental time scale (in fact they can 
even evolve on the scale of geological eras). The high degree of complexity of such systems 
makes the description of their dynamics really awkward, since apparently all the history 
of a sample has to be known in order to predict its evolution. Conversely, a useful theory 
should be able to provide a description of the behavior of the system in terms of few and 
essential effective parameters. The effective temperature, defined on the basis of fluctuation- 
dissipation relations between correlation and response functions [1], has been proved very 
fruitful in this direction (at least for mean-field models) and it is currently under intensive 
experimental investigation [2]. 

To introduce the concept let us consider the following experiment. Prepare a system (e.g., 
a glass, a ferromagnet, etc.) in an equilibrium state corresponding to a high temperature 
T (where "high" means, here, greater than any critical or glass transition temperature). 
At time t = 0, quench the system to some temperature T < T by taking it into contact 
with a thermal bath at temperature T, and let it evolve. On a general basis, one expects 
that the relaxation towards the equilibrium state corresponding to T is characterized by 
two different regimes: (A) a transient one with non-equilibrium evolution, for t < t eq (T), 
and (B) a stationary regime with equilibrium evolution for t > t cq (T), where t eq (T) is some 
characteristic equilibration time of the system. During (A) the behavior of the system 
is expected to depend on the specific initial conditions and both time-reversal symmetry 
and time-translation invariance are broken, while they are recovered in regime (B): The 
dynamics of fluctuations is given by the "equilibrium" one. We are concerned here with 
those systems for which the regime (B) is never achieved during experimental times, i.e., 
for all practical purposes, t eq = oo. In this case, the standard concepts of equilibrium 
statistical mechanics do not apply and in particular two-time quantities, such as the response 
Ro(t,s) = ( gh (K ) and the correlation functions Co(t,s) = (0(t)0(s)) (O is some 

observable and ho its conjugated field, we will be more explicit later on) depend separately 
on s (usually called the "age" of the system, being the time spent in the phase with t C q = oo) 
and t, even for long times. This behavior is usually referred to as aging [3,4]. Such a useful 
quantity as the temperature T of the system is not defined in this genuine non-equilibrium 
regime. On the other hand one can address the question whether some effective temperature 
T e ff (in general different from T) can be still defined and used to understand the physics of 
the system. 

In equilibrium [regime (B)], correlation and response functions depend only on time 
differences and the fluctuation-dissipation theorem (FDT) states that 

TRo(t-s) = d 8 C (t-s), (1) 



[here and in the following we assume t > s given that causality implies Ro(t < s,s) — 0] 
where T is expressed in units ks — 1. Whatever the regime is, one can always define the 
so-called fluctuation-dissipation ratio (FDR) as [5] 

_ TR {t,s) 

Xo{t ' s) -dMty)- (2) 

As a consequence of Eq. (1), Xo(t, s) = 1 whenever £ > s ^> t cq (T). Some sort of effective 
temperature [1] in the aging regime [t eq (T) = oo] can be defined via the long-time limit of 
X c (t,s) 

Xg = lim ]im X (t,s) , (3) 

through the relation T e g = T/Xq (that reduces to the thermodynamic temperature of the 
thermal bath when equilibrium is asymptotically reached). Obviously this definition has to 
be regarded as formal as long as one is not able to establish a link between T e g and some 
thermodynamic properties. Nevertheless, in Ref. [1] it has been shown that T e g plays the 
same role as the thermodynamic temperature, in the sense that it determines the direction 
of heat flows and acts as a criterion for thermalization. Moreover, a thermometer coupled 
to the observable O measures (on a proper time scale) the temperature T e g [1]. We also 
mention that it has been argued that Xo(t, s) establishes a bridge between the dynamically 
inaccessible equilibrium state and the asymptotic dynamics for large times [6]. (See also 
Refs. [7,8] for a discussion of other properties of Xo(t, s) and T c g.) 

It has been stressed several times in the literature (see, e.g., Ref. [9]) that the effective 
temperature can be of interest in order to devise some thermodynamics for the system 
provided that its value is independent of the observable used to define it. This has been 
explicitly verified for infinite-range (mean-field) glass models [1]. Beyond these cases, the 
observable dependence of T c g has been investigated analytically for the trap model [10], 
for the one-dimensional Ising model [11,12], and for the d-dimensional spherical model [13], 
whereas numerical studies have addressed the problem for supercooled liquids [14] and for 
the two-dimensional Ising model [11]. 

In this paper we discuss the problem of the observable independence of the FDR (and 
consequently of T c g) in a simpler (compared to glasses) class of slow-relaxing systems: Crit- 
ical systems quenched from a high-temperature phase to the critical point and evolving 
according to a purely dissipative dynamics. In fact, soon after the introduction of the FDR 
Eq. (2), it was pointed out [15] that also these systems display slow-relaxation (due to 
t eq ~ £ z , where £ is the correlation length, diverging at the critical point, and z the dynami- 
cal critical exponent [16]) and aging. The FDR of the order parameter was then determined 
for a Gaussian model and for the random walk. Subsequent analytical and numerical calcu- 
lations on realistic models confirmed and generalized this picture [17-42,11-13] (for a review, 
see Ref. [43]). For our analysis we can take advantage of the powerful tools of Renormal- 
ization Group (RG) and field theory to provide analytical predictions for some FDR's in an 
e-expansion, where e = 4 — d and d is the spatial dimensionality of the system. This kind 
of study allows us to clarify once and for all whether a unique effective temperature can be 
defined for such systems in the long-time limit. 



The outline of the paper is the following. In Sec. II we introduce the model and the field- 
theoretical approach to non-equilibrium dynamics. We derive, by means of RG equations, 
the scaling forms of two quadratic (in the order parameter) observables. In Sec. Ill we argue 
that the FDR Xq (and thus T e g) does not depend on the specific observable O within the 
Gaussian (mean-field) approximation. The effects of the interaction are taken into account 
in Sec. IV, where we show by an explicit calculation that the FDR's of two quadratic 
observables differ already at the lowest order beyond the Gaussian approximation. In Sec. 
V our results are carefully compared with the numerical and analytic ones available in the 
literature. In Sec. VI we summarize the results obtained and their implications. In the 
Appendix we report the details of the computations of the Feynman diagrams. 

II. THE MODEL 

One among the simplest non-trivial models displaying slow-relaxation and aging is a 
lattice spin model in d dimensions with O(N) symmetry and short-range interactions evolv- 
ing according to a purely dissipative dynamics after a quench to the critical point. In the 
simplest instance its lattice Hamiltonian is given by 

W = -5>. Sj , (4) 

where Si is a iV-component spin located at the lattice site i, with Si 2 = 1. The sum runs 
over all pairs (ij) of nearest-neighbor lattice sites. A purely dissipative dynamics for this 
model proceeds by elementary moves that amount to random changes in the direction of 
the spin S{ (spin-flip sampling). The transition rates can be arbitrarily chosen provided that 
the detailed-balance condition is satisfied. For analytical studies the most suited are the 
Glauber ones [44], which allow exact solutions in the one-dimensional case [17-19,11,12]. 

Despite their simplicity, these models are not exactly solvable (for arbitrary N) in phys- 
ical dimensions d = 2,3, and to obtain information about the non-equilibrium critical dy- 
namics, one has to resort to numerical simulations [20,21,11,29-33]. 

To investigate analytically the dynamical behavior in physical dimensions, we take ad- 
vantage of the universality considering the time evolution of a iV-component field y?(x, t) with 
a purely dissipative dynamics (Model A of Ref. [16]). This is described by the stochastic 
Langevin equation 

^(x,t) = -fi-^M- + 6(x,t), (5) 

6<pi(x,t) 

where Q is the kinetic coefficient, £(x, t) a zero-mean stochastic Gaussian noise with 

(&(x,t)&(x',0> =2Q6(x-x!)6(t-1?)6 ij , (6) 

and Ti.[tp] is the static Hamiltonian. It may be assumed, near the critical point, of the 
Landau-Ginzburg form 



HM= / d d x 



_(V^) 2 + -7 W > 2 + -<?0</ 



(7) 



where r oc T is the temperature parameter, assuming its critical value r 0iC for T = T c 
( r o,c — 0, within the analytical approach discussed below), and Qq is the bare coupling 
constant of the theory. This coarse-grained continuum dynamics is expected to be in the 
same universality class as the lattice models with O(N) symmetry, short-range interactions, 
and spin-flip dynamics [16]. 

The equilibrium correlation and response functions can be obtained by means of the 
field-theoretical action [45,46] 



S[(p,<p] 



dt / d x 



'.dtp SH[p] 

<p m +n<p ^ p - 



ifVtip 



(8) 



Here <p(x.,t) is a iV-component auxiliary field, conjugated to the external field h in such a 
way that 7i[tp, h] = Ti.[tp] — f d d xh<p. As a consequence, the linear response to the field h of 
a generic observable O is given by 



S(Q) 
6hi(x,s) 



Q{<Pi(-x, s)0) , i 



N . 



0) 



For this reason <p(x.,t) is termed response field. 

Within this field-theoretical formalism it is possible to show that the FDT holds for 
generic Rq and Co. This has an illuminating derivation in a supersymmetric formulation, 
where the FDT's are the Ward identities due to the supersymmetry [45]. 

The effect of a macroscopic initial condition y?o( x ) — y?( x , ^ — 0) may be taken into 
account by averaging over the initial configuration with a weight e~ H °^ where [47] 



H [(po 



^^[v9 (x)-a(x)] 2 , 



(10) 



that specifies an initial state o(x) with Gaussian short-range correlations proportional to 

Following standard methods [45,46] the response and correlation functions may be ob- 
tained by a perturbative expansion of the functional weight e~^ '■ (p '^' +H °^' Po '' in terms of the 
coupling constant go. The propagators (Gaussian two-point functions of the fields (p and (p) 
of the resulting theory are [47] 



(<Pi{q, s)<pj(-q, t)) = SijRq(t, s) = 5 tj 0(t - s)G(t - s) 



{<Pi(q, s)<pj(-q, t)) = $ijC°(t, s) 



Sij 



<r + r 



G(\t-s\) + 



r + q_ 2 



To 



l)G(t + s) 



where 6(t) is the step function [6(t < 0) = 0, 9{t > 0) = 1] and 



G(t) 



-fl(^ 2 +r )t 



(11) 

(12) 



(13) 



The response function Eq. (11) is the same as in equilibrium. Eq. (12), instead, reduces to 
the equilibrium form when q^ and both times t and s go to infinity while r = t — s is kept 
fixed. In the following we will assume the Ito prescription (see, e.g., Refs. [48,45]) to deal 
with the ambiguities arising in formal manipulations of stochastic equations. Consequently, 



all the diagrams with loops of response propagators have to be omitted. This ensures that 
causality holds in the perturbative expansion [47,48,46]. From the technical point of view, 
the breaking of time-translation invariance does not allow the factorization of connected 
correlation functions in terms of one-particle irreducible ones as usually done when time- 
translation invariance holds. As a consequence, as it is the case when dealing with surface 
critical phenomena [49], all the computation has to be done in terms of connected functions 
only [47]. Furthermore it has been shown [47] that Tq 1 is an irrelevant variable for the RG 
flow affecting only the correction to the leading long-time scaling behavior we are interested 
in. In view of that we fix it to its fixed-point value Tq 1 = from the very beginning of the 
calculation. 

From scaling arguments [20], and more rigorously from the solution of RG equations [47], 
it is known that the zero-momentum response and correlation functions of the basic fields 
satisfy the scaling forms (see Ref. [43] for a review): 

i? q=0 (t, s) =A R (t- s) a (t/s) e F R {s/t) , (14) 

C q=0 (£, s) = A c s(t - s) a (t/s) e F c ( S /t) , (15) 

where a = (2 — r] — z)/z, z is the dynamical critical exponent, rj the anomalous dimension 
of the fields [45], and 9 the initial-slip exponent [47,48]. We single out explicitly the non- 
universal amplitudes A Ry c by fixing -Fr,c(0) = 1. With this normalization F Ry c are universal 
scaling functions. From the previous scaling forms one deduces that 

<9 s C q=0 (i, s) = A dc (t - s) a (t/s) e F dc (s/t) , (16) 

where the non-universal amplitude Age has been defined so that Fqc{0) = 1. Accordingly 
one has A dc = A c (l — 9). 

Using Eqs. (14) and (16) one finds that 

x ( . x- gg=o(M) A R F R (s/t) 

^q=0{t, S = — — — — — , [if 

d s C q= o{t,s) A dc F dc {s/t) 

is a universal amplitude-ratio (in the sense of Ref. [50]) being the ratio of two quantities 
[Rq = o(t,s) and <9 s C q= o(£, s)} that have the same scaling dimensions. Furthermore it is a 
function of the ratio s/t only, and not of s and t separately 1 In lattice simulations (especially 
in the literature concerning glassy systems [7,8]) response and correlation functions are often 
measured in the real space x, instead of in the momentum space q as done here. From the 
scaling forms given above one can derive the analogous ones in the real space for R x , C x , 
and then define X x= o = R x =o/d s C x= o (originally introduced in Refs. [5,15]) in analogy with 
Aq =0 . In general one expects A x=0 ^ Aq =0 . Nonetheless it has been argued [25] that the 
long-time limit of the universal amplitude ratio 

X™ = lim lim /^ =o(t ; S) N = lim X q - (s/t) = 4^ = , f* ^ , (18) 

m ^oo^oo0 a C q=o (M) -A-o q " 0W ; A dc Ac0--0) 



1 This is an important difference compared to mean-field glassy model, where, instead, it turns 
out that X x= o in the long-time regime can be written as a function of C(t, s) (see, e.g., Refs. [7,8]). 
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(the subscript M refers to the fact that y?q=o oc M, the average magnetization) is equal 
to the same limit of X x=0 . This equality was also confirmed by numerical simulations [11] 
(even if the numerics of Ref. [11] have been questioned [51]). 

A. Scaling forms of composite operators 

We now derive scaling forms analogous to Eqs. (14) and (15) for the correlation and 
response functions of local composite operators, focusing on those of the form (p m . However, 
the derivation is completely general and can be easily applied to any other operators. 

Let us consider an observable O (a composite operator, using the field-theoretical termi- 
nology) having h as a conjugate field (e.g., O is the energy density and ho the temperature) 
and coupling to H according to TL i— ► 7i + h 0. As a consequence, the dynamical functional 
S changes according to S \— ► So — S + QhoO, where the associated operator O is given by 

6 = /**.**,.)^. (19) 

The linear response of an observable A to a variation in the field ho can be expressed as 

HA) ho 



Sh 



o 



= n{AO) , (20) 

h o =0 



where {-)h stands for the average over the dynamics associated with the dynamical func- 
tional in the presence of ho- This generalizes Eq. (9). 

To render finite the correlation functions with insertions of the operator O, of the form 

{[^nmorpf) , (21) 

one additional renormalization (compared to those necessary when O and O are not inserted, 
evaluated in Ref. [46]) is required: Ob = Z Or. Here and in the following with the subscript 
B we indicate the bare quantities and with R the renormalized ones, whose correlation 
functions are finite upon removing the regularization [45]. In the case of operators mixing 
under renormalization, O has to be understood as a vector of suitable operators, while Zo 
will be in general a renormalization matrix [45]. The presence of additive renormalizations 
does not change the scaling arguments presented below. In view of Eq. (19) one finds that 
& B = Z d O R with 

Z d = (Z^/Z^Zo , (22) 

where Z 9 and Z^ are the renormalization constants of the fields, defined, as usual, by <ps — 
Zy 2 Lpn and (pB = Z-, <pr. The correlation function (21) can be renormalized according to 

M^UO]°PT)b = z;/ 2 zy 2 z z d d ([ v r[mo]°[df) R , (23) 

where, on the r.h.s., (•)# means that all the bare quantities have been replaced by the 
corresponding renormalized ones. Applying standard techniques it is possible to write the 



RG equations by introducing appropriate RG functions [45]. For the theory defined by So, 
in addition to the RG functions of the theory with action S, one has to introduce two new 
functions go = f^d^ln Zo\o and q-q = fid^lnZ^lo 2 (fi is the scale at which the theory has 
been renormalized [45]). With |o we indicate that the differentiation has to be done with 
fixed bare parameters. In view of Eq. (22), go and g^ are related by qq = go + (g<p — £v)/2 
(where g v and g^ are the usual RG functions for the fields). Combining dimensional analysis 
with the solution of the RG equations one finds that the scaling dimension of the correlation 
function ([(9]°[C , ]°) R in the (x, £)-space is given, at the infrared fixed point of the theory, by 

5(o,5) = o[0} scal + 5[d} scal . (24) 

With [-] sca i we indicate the scaling dimensions. In terms of the fixed-point values of the 
RG functions and the canonical (engineering) mass dimensions of the operators (denoted by 
[•]can) they are expressed as: 

[fiscal = [0} cw + Q0(g*), 

[0] scal = [<9] ca „ + g d (g*) , (25) 

where g* is the fixed-point value of the renormalized coupling constant. In analogy with 
the anomalous dimensions i] = g v (g*) and fj = g^(g*) (such that [<p] sca \ — (d — 2 + rj)/2 and 
[fiscal — (d + 2 + fj)/2 [46]) one introduces rjo and rj^: 

d -±^ - [0],„ . (26) 

(Note that, in contrast to rj and fj, r\o and rj^ do not generally vanish in the free theory 
g* = 0.) It is easy to verify that, in terms of r]o, the critical exponent 70 of the susceptibility 
((£>£>) ~ \T - T c \~i°) is given by 

1o = i/(2 - r)o) , (27) 

where v is the critical exponent of the correlation length. As a consequence of Eq. (19) 
one has [0] can = [C] can + [ip] can - [ip] cmi = [C] can + 2 (recall that in Model A dynamics 
Mean — (d — 2)/2 and [<p] can — (d + 2)/2). Using (fj — f])j2 = z — 2 (a consequence of the 
fluctuation-dissipation theorem [45]), we have 

V °- Vo = Z -2 (28) 

and [0} scal = [C] sca i + z, leading to 



2 go and q~ are the anomalous dimensions of O and O, respectively. In the literature they are 
usually referred to as 770 or 70 [45]. However, in the present case we use r\o and 70 to indicate 
critical exponents. 



5(2,0) = d-2 + Vo , 

6(1,1) = d-2 + rio + z. (29) 

Taking into account that [time] sca i = —z, one finds (hereafter we set Q = 1) 

(0(q, t)0(-q, s)) = (t- s) a ° +1 F c (q z (t - s), s/t) , (30) 

(0(q, t)0(-q, s)) = (t- s) a °F R (q z (t - s), s/t) , (31) 

where 

_ 6(1,1) -d 2-rio-z 
ao = = • \6Z) 

z z 

Let us now discuss the effect of the temporal surface. The scaling dimension of (po has 
been computed in Ref. [47], i.e., [y5 ] S cai = [y5] SC ai + ^o/2. (In terms of i] Q and z the initial- 
slip exponent is given by 9 = —f]o/(2z) [47,48].) Consider now the case of the observables 
0^ m \t) constructed with m fields ip at time t, i.e. 0^ m \t) ~ <p m (t). Eq. (19) implies that 
(m) (t) ~ ^(t)v? m ~ 1 (t). The observable (m) , obtained from £> (m) by replacing all the fields 
(p with those tp, has dimension given by 

[O^U^lO^U^ + mz. (33) 

Moreover, keeping in mind the scaling dimension of (po one has, for t — 0, 

[^ m) ]scai=[0^] scal -m^, (34) 

where O^ = (m) (* = 0). 

When inserted into correlation functions, (dt<p)t=o = 2Q(po, while the insertion of an 
initial field ipo vanishes [49]. Thus, for the operators CK m ) and O^ we expect the following 
short-distance expansion for t — > 

C (m) (t) ~ p(t)Cf% + h.o.c.f. , (35) 

a (m) (t)~p(t)^™ + h.o.c.f., (36) 

where h.o.c.f. stands for higher order composite fields that contribute to this expansion only 
with subleading terms. As a consequence, the scaling dimensions of p and p are given by 

[p] S cai = -mz + m9z and [p] SC£d = -(m - l)z + m6z, (37) 

where we used Eqs. (33), (34), and the FDT. Taking into account that [time] sca i = —z, one 
concludes that, for s — ► 0, 

(m) (s) ~ s m - me and (m) (s) ~ s m -^ me . (38) 

It is now possible to rewrite the scaling forms (30) and (31) in a way that shows explicitly 
the behavior of the scaling functions Fq and Fr for s — ► 0, i.e., 

(0 (m) (q,t)0 (m) (-q,s)) = s(t - s) ao (t/s)-( m -V +me Fg(q z (t - s),s/t) , (39) 

(0 (m) (q,t)d (m) (-q,s)) = (t-s) ao (t/s)-^ m -^ +me Fg(q z (t-s),s/t), (40) 



where now the functions F® and F R are regular for s — > 0. Furthermore they are also 
universal once we fix the normalization for small arguments. 

Let us focus on the scaling properties of the correlation and response functions of the 
observables with m — 2. Because of the O(N) symmetry of the underling theory only two 
quadratic zero- momentum operators with different scaling dimensions exist, that can be 
written as [45] 

- d d q N 



E(t) ee [d d x £ y^(x, t) = (-*- £ <ft(q, *)<ft(-q, t) , (41) 

J i=i J \ zn ) i=i 

T^it) = jd d x ^(x,t)^(x,t) - ^SijElt) = j-^- <Pi(q,t)<pj(-q,t) - ^E(t) ■ ( 42 ) 



The corresponding response operators are given by 
N '■ d d q 



E (t) = E J?A- d 2 ^(q, t)<pi(-q, t) , (43) 



i=i- (2^) 



Ty(t) = /(^7 [^(q^)^-(-q,i) + ^(q,«)^-(-q,«)] - ^<W) • ( 44 ) 

In the following we will refer to them, generically, as 0(t). The corresponding scaling 
functions are 

C°(t, S ) = (0(t)0(s))=A° S (t- S ) a v(t/ S )- 1+2e F$( S /t) , (45) 

R°(t,s) ee (0(t)d( 8 ))=A°(t-8) a °(t/s)- 1+ » f F°(8/t) , (46) 

where the non- universal amplitudes A® R are determined so that Fq r (0) = 1. 
In terms of these quantities we can write the (zero-momentum) FDR as 

*°fc" s S& (47) 

(here we restore the actual value of fi, currently set to 1) that, as its analogous (17), depends 
only on the ratio s/t and it is a universal function. In particular its universal long-time limit 
reads 

1 A° 

^ = &&*><*■«> = 2(1^- < 48 > 



1. Overview of the known values of the exponents 

For the two specific cases of the quadratic operators E(t) and Tij(t), the exponent ao 
can be expressed in terms of more familiar critical exponents. In fact we obtained ao = 
(2 — i]o— z)/z, with rjo scaling dimension of the operator O. Using Eq. (27) we can express 
it in terms of the susceptibility exponent, given by 3 (see, e.g., Ref. [52]) 



3 When hyperscaling holds (i.e., for d < 4) 2 — dv can be replaced by a, the specific heat exponent. 
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loA 2 - dV ' ° = E ' (49) 

where 0t is the so-called quadratic crossover exponent [45,52]. The exponents 7^ and 7t 
are exactly known in d = 2 [45]. Even in d = 3 they are known with very high accuracy, in 
fact they have been computed up to five loops in the e-expansion [53] and six (sometimes 
seven) loops in fixed dimension d = 3 [54,52] for generic values of N (see Ref. [55] for a 
review) . The dynamic critical exponents are currently known with a much less accuracy. In 
fact z has been computed up to three loops in the e-expansion [56] and up to four in fixed 
dimensions [57], whereas in known only up to two loops in the e-expansion [47]. For d > 4 
the mean-field results hold: ao = 1 — djl (for both E and T) and 8 = 0. 

For later convenience, we report here explicitly the 0(e) expansion of a^ 4 (recall that 
z = 2 + 0(e 2 ) [16]) 



a a 



4 — N 
■1 + -— ^ e + 0(e 2 ), = E 



2(N + 
2(N + 



-l + ^ + l e + 0(e 2 ), O 



(50) 



and 9 

The exponents are exactly known in the limit iV — > 00: z/ -1 = d — 2, -2 = 2, = 1 — d/A 
[47], and 0r = 2z/ in d < 4, whereas 1/ = 1/2, z = 2, 6 = 0, and 0^ = 1 in c? > 4 (mean-field 
exponents). Using these values, the scaling forms for the energy in d < 4 are 

C E (t, s) = A E C s(t - sY^-^t/sf-^F^s/t) , (52) 

R E (t, s) =A E R (t- s)^-\t/sf^F E {s/t) , (53) 

whereas for T in d < 4, 

C T (t, s) = A T C s{t - sf-^it/sf-^F^is/t) , (54) 

R T (t, s) =A T R {t- sy-^it/sy-^F^s/t) . (55) 

For d > 4 the scaling forms for E and T are the same: 

C°(t, s) = A% s{t - sy- d / 2 (s/t)F°(s/t) , (56) 

R°(t, s) =A°(t- sy- d l\s/t)Fi{s/t) . (57) 



4 In passing, let us mention that, for N = 1, (fa has a non-trivial value, even if the operator Tjj(t) 
is not defined for the Ising model. This fact has an interpretation in terms of a gas of N— color 
interacting loops (see, e.g., Ref. [58]) belonging to the same universality class as the O(N) model. 
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III. THE GAUSSIAN APPROXIMATION 

For the Gaussian model the response and correlation functions are known exactly, so we 
can evaluate the FDR (in Ref. [15] this has been done directly in real space). From Eqs. (11), 
(12), and the definition (17), properly generalized to q 7^ 0, one finds [25] 

Q/?° 1 

X ^ t ' ^ = d^C® = 1 + e -2n(q2+r )s • ( 58 ) 

If the theory is non-critical (tq 7^ 0) the limit of this ratio for s — > 00 is 1 for all the values of 
q, in agreement with the idea that in the high-temperature phase all the fluctuating modes 
have a finite equilibration time, so that equilibrium is recovered and the FDT applies. In the 
critical theory the limit ratio is again equal to one when q^ 0, whereas for q = one has 
X® =0 (t,s) = 1/2. This shows that the only mode that "does not relax" to the equilibrium 
is the zero mode in the critical limit. This picture (already presented in Ref. [25]) has been 
confirmed by a one- loop computation [25]. 

In the case of the Gaussian theory, with go = in the Hamiltonian (7), the FDR can be 
easily computed for a generic observable. In particular we show that the FDR for a set of 
local one-point observables is always equal to 1/2 in the long-time limit. This allows for a 
definition of a unique effective temperature in the Gaussian model. We consider operators 
of the form 0j„ = d l (p n . The correlation and response functions of all local operators can 
be written in terms of those of Oi jU . 

The critical (i.e., r$ = 0) response and correlation functions of the order parameter 
Eqs. (11) and (12) are given, in real space, by (their diagrammatic representation is reported 
in Fig. 1 (b) and (a), respectively) 

R x (t,s) = F d (t-s,x), (59) 

C x (t, s) = K d (t - s, x) - K d (t + s, x) , (60) 

where (we set Q — 1) F d (r, x) = J(dq)e~ cl r e~* q ' x and K d (r, x) = /(dg)q~ 2 e~ q r e^ q ' x (with 
(dq) = d d q/(2ir) d ). 

Let us explain our argument considering first the operators Oq^. The two-point correla- 
tion functions of (9o,n is given by the (n — 1)— loop diagram with the two points connected 
by n correlation lines, as shown in Fig. 1 (c). 

In the real space its expression is simply given by the product of n correlators. Thus 

C?(t, S ) = c n [C x (t,s)] n , (61) 

(c n is the combinatorical factor associated with the diagram) whose derivative is 

d 8 C?(t, s) = c n n[C x (t, s)\ n - l d s C x (ti s) . (62) 

Analogously the response function is given by the diagram depicted in Fig. 1 (d), obtained 
from that one contributing to the correlation function (Fig. 1 (c)) by replacing an order- 
parameter correlator with a response function: 

R°(t, s) = c n n[C x (t, s)} n ^R x (t, s) , (63) 
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C x (t,s) 
(s,0) (t,x) 

(a) 



Rx{t,s) 

(s,0) > (t,x) 

(b) 




n 



(i,x) (s,0) 




n 



(c) (d) 

FIG. 1. Diagrammatic elements: (a) correlation and (b) response functions of the or- 
der parameter. (c) Diagrammatic representation for the two-point correlation function 
C®(t,s) = (Oo >n (t, x)C\), n (s,0)} and (d) response function R®(t,s) = (Oo jn (t,x.)Oo, n {s,0)}. See 
the text for further explanations. 



where the factor n comes from the fact that the response operator is Oo, n — nOo,n-i<P- Note 
that the combinatorial factor c n is the same as for the correlation function. 
The FDR is given by 



X^(t,s) 



i?x°'"(t,s) 



d s C°°* n (t,s) d,C x (t,s) 



X?(t,8) 



(64) 



Thus we obtain the remarkable result that the FDR of powers of the field in real space is 
equal, for all times, to the FDR of the field. 

Before considering the effects of the derivatives let us consider the previous relation in 
momentum space. Remembering that the product of two functions after a Fourier transfor- 
mation becomes a convolution, one has 



R°(t,s)=c n n(C*---*C*R) (l , 



d s C^(t,s)=c n n(C 



*C*C" 



q • 



(65) 
(66) 



where * is the convolution, C and R are in momentum space (with the time dependence 
understood), ••• means n — 1 times, and C = <9 s C q . For q = the previous relations 
become 



i2?=o(*. s ) = c n nj(dp)(C * • • • * C) P R_ P , 
d s C% =0 (t, s) = c n n f(dp)(C * • • • * C) p C'_ p . 



Thus 



X£o(t,8) 



-1 



J(dp)(C*---*C) p C'_ p J(dp)(C * • • • * C) p R_ p Xz 
J(dp)(C*---*C) p R_ p - 
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-i 
p 



J(dp)(C 



C) P R. 



(67) 
(68) 

(69) 




B c (tuh) 




Gi(t\,t2 




L 4 (h,t 2 ) 

FIG. 2. Diagrams contributing to the correlation function. Gi(t,s) and L±{t,s) have also 
symmetric counterparts that have been considered. 

i.e., the inverse of X® =Q is a weighted average of XZp, with weight (C * • • • * C) P R- P , 
that in the limit t — > oo, s — »■ oo in the proper order, is peaked around p = 0, giving 
A^ = lim s _ oo lim t _ +oo A'^ =0 (i,s) = lim^oo lim^^ X p=0 (t, s) = X%}. Note that in momen- 
tum space, at variance with the relation Eq. (64) holding between FDR in real space, only 
the long-time limit of the FDR reproduces the FDR for the fields. This is in agreement with 
some explicit calculations we made for (p 2 (see below) and y? 3 FDR's. 5 

Let us now take into account the effect of the derivatives. In momentum space they 
amount simply to a multiplication by q* and so they change the weight in Eq. (69) by a 
factor q 2i (one q* for each insertion), not affecting our conclusion on the long-time limit. 

Let us note that if instead we introduce some non-local operators, as e.g. e d ° (with O 
a local operator), the weight in Eq. (69) is no longer exponential in q 2 and the long-time 
limit of FDR may be different from 1/2. 



IV. TWO-LOOP COMPUTATION 

Here we present the details of the two- loop (i.e., up to 0(e) in the e-expansion) pertur- 
bative computation of the correlation and response functions for the zero-momentum energy 
and the tensor operators given by Eqs. (41) and (42). The one-loop results are the Gaussian 
ones obtained in the previous section. 



5 The Gaussian </? 3 FDR is obtained from Eqs. (A9) and (B6) of Ref. [26]. The calculation is 
straightforward, but cumbersome. 
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B R (t,s) 




G 2 (t,s) 





Li{t,s) 



t s< 




L 2 (t,s) 



L 3 (t,s) 



FIG. 3. Diagrams contributing to the response function. 

The diagrammatic contributions are reported in Fig. 2 for the correlation function 
C (t, s) and in Fig. 3 for the response function R°(t, s). In Table I we list the symmetry 
factors of the diagrams (depending both on the global topology and on the external legs) 
and the corresponding color factors (depending on global topology and the index structure 
of the vertices). 

For later convenience let us introduce the function 



In(t) 



d d q e 



2q 2 t 



(27r) d q 



d n 2n 



In{l)t 



-d/2+n 



where 



I n (l)=N d 2 n - d l 2 - 1 T(d/2 



n) 



(70) 



(71) 



and N d = 2/[(Aii) d l 2 T(d/2)}. 

Thus we obtain, for the one-loop diagrams 6 



6 Note that without the temporal surface at t = 0, i.e., using the equilibrium correlator C^ eq '(ti — 
£2), -Be would be given only by ^(l^i — *2 1 ) • For d = 4 — e, I2 has a pole [^(l) oc T(2 — d/2)], that 
is usually subtracted by introducing the well-known additive renormalization for the tp 2 correlation 
function. In this case, however, the time-boundary term yields two additional factors that exactly 
cancel the pole of the first term. Accordingly, in contrast to the equilibrium case, no additive 
renormalization has to be introduced to render finite the energy correlation function. Thus the 
equilibrium renormalized correlation functions for q = cannot be simply recovered by taking the 
limit t, s — > 00 with t — s fixed, as usually done for field correlation functions. 
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TABLE I. Symmetry and color factors of the Feynman diagrams depicted in Figs. 2 and 3. 



Symmetry 
Corr. Resp. 






Color 
coeff. E 


T 


Br, B C 

G\i G2 

Li,L,2, L3, L4 


2 

2 

2 


2 
2 
1 






d= N 
C 2a = N(N + 2)/3 
C 2b = 7V(AT + 2)/3 


1/2 

1/3 

(N + 2)/6 


B c (t 1 ,t 2 ) 

B R (t, s) 


" /(2vr)^ (tl 


i)R% s) - 


(l*i 
= h 


(* 


t 2 |) - 2/ 2 (max{£ 1 , t 2 }) + / 2 (ti + t 2 ) 
-s)-h(t), 


, (72) 
(73) 



where here and hereafter we set Q = 1 to lighten the notation and t > s (for t < s the 
diagrammatic contributions are zero due to the causality of the response functions). 

In terms of the two-loop Feynman integrals the response and correlation functions read 

E°{t, s) = 4dS R (i, s) - 4C 2a gG 2 (t, s) - 2C 2b g[L 1 {t, s) + L 2 (t, s) + L 3 (t, s)] + 0(g 2 ) , (74) 
Cg(t, s) = 2C7iS c (t, s) - 2C 2a g[G 1 (t, s) + G^s, t)] - 2C 2b g[L 4 (t, s) + L 4 (s, t)} + 0(g 2 ) . (75) 

The calculations of Gi and Lj are quite cumbersome and all the details are reported in the 
Appendix. Taking into account the expressions therein we find renormalized response and 
correlation functions that are in agreement with equations (46) and (45), with the exponents 
given in Eqs. (50) and (51). The non-universal amplitudes and scaling functions are given 
by (as usual we introduce g = N^g to simplify the notation) 



A 



E,T 

R 



, e . /e C 2a g\ ( . , x . C 2b g 



N d C x 



F R (X) 



A 



E,T 

C 



F C ( x ) 



e 

1 -- 

2 

1 — x 



C 2 b 9 

d 4 

e 



lofffj 



X 



logfl — x) 



X 



C 2a g\ , . ■. C 2 bg C 2a 5 _' 
2-d2) {lE + Y ° g2)+ d2 + d2-A 9 . 



N d d 



x- 



-log(l-x 2 )--Wi(x) 



C*2a J,,,, n C 2 b g TX/ 1 x 
dl W ^ x) - d~2 W3{x) 



(76) 

(77) 
(78) 
(79) 



where the functions Wi(x) have been defined in Eqs. (A24), (A25), and (A26). The previous 
equations are valid up to 0(g 2 ,eg,e 2 ). 

According to Eq. (48), Xq can be expressed as an amplitude ratio 



voo 



K 



2Ag(l - 9) ■ 
Using the non-universal constants given in Eqs. (76) and (78) we find that 



(80) 






1- 



C26 g* C 2a 5 
~d~8 ~ ~d24 



9 



0(e 2 



(81) 
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where 

v" = N^S € + 0{t2) - (82) 

is the fixed-point value of the coupling constant [45] . Taking into account the combinatorical 
factors in Table I we finally find: 






2 V 3A + 8 



I 2 V 12 N + 8 ' 



(83) 



Recall that the one-loop FDR for the order parameter (magnetization) is [25] 

From these results we conclude that the long-time limit of the fluctuation-dissipation ratio 
depends on the particular observable chosen to compute it. 

As a final remark let us comment on the connection between the correlation functions 
of if 2 in Model A dynamics and those of the conserved density e in the associated Model C 
[16]. In equilibrium dynamics it is usually easier to compute (y?V 2 ) i n terms of (ee) (see, 
e.g., Ref. [59]). However, this is no longer possible when studying non-equilibrium dynamics 
(considered in Ref. [28]) because of the connection between (p and e does not carry over to 
this case. 



V. COMPARISON WITH OTHER RESULTS 

A. The spherical model and the N — > oo limit 

The static critical behavior of the O(N) model in the limit N — *■ oo is known to be 
equivalent to that of the spherical model, defined by the Hamiltonian 

where s^ are real numbers subjected to the constraint Xa s f = L d (L being the linear di- 
mension of the d- dimensional lattice, assumed for simplicity to be hypercubic), and the sum 
runs over all nearest-neighbor pairs (ij). Indeed Stanley proved [60] that the free energies of 
the two models are exactly the same. From this equality it follows that critical exponents, 
scaling functions etc. are equal. The same equivalence also holds for equilibrium critical 
dynamics defined in the spherical model by means of a Langevin equation as Eq. (5) 7 . As 



7 This equation has to be properly modified in order to prescribe a dynamics that is compatible 
with the spherical constraint. 
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far as we are aware, this equivalence has not yet been carried over to the non-equilibrium 
critical dynamics we are interested in. Let us recall that the free energies of the models in 
the presence of a spatial boundary (whose corresponding field theory looks very similar to 
the non-equilibrium one that we are considering) are well known to be not equal [61]. 

The spherical model has attracted a lot of attentions, since its essential Gaussian Hamil- 
tonian makes it exactly solvable although the resulting critical behavior is not mean-field 
like because of the spherical constraint. The FDR of the magnetization was calculated 
by Godreche and Luck, who obtained X^ = 1 — 2/d. This result is compatible with the 
two-loop e-expansion [26] for N = oo and, moreover, it agrees with the exact result we are 
going to derive. Recently Sollich [13] has determined several FDR's of quadratic (in the 
spin Si) operators that could be compared with our results. He considered the bond energy 
observable B[ = |(si — Sj) 2 , the product observable Pi = S[Sj (with i, j nearest neighbors) 
both in the real and in the momentum space, and the total energy. The exact results of his 
analysis are that Xp = Xpf = X^ ^ X|? (the expansion of X|? close to four dimensions for 
d < 4 nicely agrees with our two- loop e-expansion for iV = oo). The scaling forms for these 
observables have been also derived and they agree with our equations (31) and (30) with 
as = — 1 — d/2, dp = 1 — d/2, and ap = d/2 — 3 in d < 4 (and ap = ap = ap + 2 = 1 — d/2 in 
d > 4). All these findings agree with our calculation apart from the fact that the observable 
P of the spherical model cannot be naively identified with (p 2 for N = oo. On the other 
hand ap = ap — 2 agrees with the naive identification of B as the Laplacian of P in the 
continuum limit. 

We report now the calculation of some FDR's for the O(N) model directly for N = oo, 
using the well-known property that for N = oo the fourth order interaction term can be self- 
consistently decoupled g/N({p 2 ) 2 — > gC^=o(t)ip 2 (here go is replaced by g/N) [45]. Taking 
advantage of this decoupling the exact response function of the theory (for 2 < d < 4) has 
been computed [47], finding: 

RM,s)=9(t-s)(-j e- q ^, (86) 

whereas, for the correlation function [47] 

CJt, s) = 2 dt'RJt, t')RJs, t') = 2(ts) 1 - d/4 e-* 2(t+s ^ / dt' ^l^e 2 ^' . (87) 

Jo Jo 

In particular for q = the previous expressions become 

/t\ 1 ~ d/4 2 /t\ 1 ~ d/4 

R q=0 (t, s) = 9(t - s) (-J and C q=0 (t, s) = d/2 _ 1 * [-) ■ (88) 

According to Eq. (18) it is straighfroward to compute the FDR for the order parameter, 
finding X^ = 1 — 2/d, the same result as for the spherical model. 

Now we consider, in the limit N — > oo, the scaling functions (45) and (46) of the response 
and the correlation of the quadratic operators. The explicit computation for 2 < d < 4 of 
the scaling functions, the non-universal amplitudes, and, eventually, the FDR of composite 
operators is not as straightforward as one could have erroneously expected. In fact, if one 
roughly assumes, on the sole basis of the decoupling in the large- N limit, that the theory 



(8,0) 





fc-1 



A: 
(a) 



fc+1 




(t,x) 



( S ,0) 





(t,x) 



(b) 
FIG. 4. Chains of bubble diagrams for the (a) correlation and (b) response function of quadratic 

observables. For the correlation (a), the index k runs between 1 and n. 



is essentially Gaussian with renormalized two-point functions given by Eqs. (86) and (87), 
the result Xq = X^ would follow from the argument we gave for the Gaussian Model, that 
makes no use of the specific expressions of the response and correlation functions. 

Let us consider more closely the case of the energy E. For the two-point function 
(E(t)E(s)), a family of diagrams (chains of bubble diagrams, see Fig. 4) that are of order 
0(N°) even in the limit N — > oo exists (although the coupling constant is of order 1/N, each 
bubble carries a combinatorical factor N). These diagrams are not accounted for by a simple 
renormalization of the two-point functions of the order parameter. This fact is well known 
for static observables (see, e.g., the calculation of the structure factors {{p 2 (q){p 2 (— q)) up to 
four loops in Ref. [52]). If and only if the contribution of such diagrams to X°° vanishes (or 
is cancelled out) one has X^ = Xjfi- The results of the previous sections indicate that this 
is not the case. 

Note that if one considers instead of E the operator T then the chains of bubble diagrams 
are depressed by a combinatorical factor of order 1/N and so they do not contribute to the 
correlation and response functions for N = oo. On this sole basis one concludes that in the 
limit N = oo, X£? = X^ to all order in e and not only at the first one as we have explicitly 
obtained (compare Eq. (83) with Eq. (84) considering the limit N — > oo). 

This analysis indicates that the observable P in the spherical model has a scaling behavior 
(and X°°) that is the same as that of T, contrarily to the naive expectation, suggesting 
instead E. 

Finally let us mention that for all the quantities we considered here, the non-equilibrium 
dynamics of the spherical and O(oo) models are exactly the same. This fact calls for a more 
rigorous investigation of a possible correspondence between the two models beyond the case 
of equilibrium dynamics. 
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FIG. 5. Xg for O = M,E,T and N = 1,2,3, and oo. For X$ we report the two-loop result 
[26]. For X|? we report the direct estimate and those obtained by constraining at the lower critical 
dimensions (with linear constraint for N = 1, 2, 3, oo and with quadratic one for N = 2, 3, oo). For 
Xj? only the direct estimate is reported. For N = oo we also report as "E exact" the exact result 
for the spherical model, from Ref. [13]. 

B. The Ising model 

Apart from the spherical model, the FDR of composite operators has been considered so 
far only for the one- and two-dimensional Ising models, analytically [11,12] and numerically 

[11]- 

In one dimension the several FDR's considered [11,12] turned out to equal X^ = 1/2, 

apart from that of the energy: X|? = 0. This fact has been interpreted in Ref. [11] as 

an interplay between criticality and coarsening, a peculiarity of those models (such as the 

one-dimensional Ising model) with T c = 0. Instead, our result indicates that XJ? < X^ is 

a more general property, rigorously true close to four dimensions (e<l) and perhaps valid 

up to e = 3. 

In two dimensions, the numerical results Xj% = 0.340(5) ~ X^ = 0.33(2) apparently 
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indicate the equality of the two FDR's [11]. Obviously our result calls for a more precise 
determination of X|? = 0.33(2) to understand whether this apparent equality is due to the 
relatively low precision of such a measure or to the fact that d = 2 is a peculiar case for 
some still unknown reasons. In any case, let us stress that X|? from the 0(e) result for 
e = 2 gives a value much smaller than X^. The direct estimate (unconstrained [u]) from 
the two-loop series in Eq. (83) for e = 2 (giving -X"|?[u] ~ 0.28) is probably unreliable as it 
was the case for Xfi. (X^[lloop] ~ 0.42 [25] whereas X^[21oop] = 0.30(5) [26].) To obtain 
a more reliable result (without computing the 0(e 2 ) term, that seems to be a very difficult 
calculation, requiring the evaluation of three-loop diagrams) one can constrain linearly [1] 
the O(e) result to assume the exactly known value for d — 1 (i.e., e = 3), as usually done for 
this kind of expression (see, e.g., Ref. [55]). Assuming a smooth behavior in e up to e = 3, 
one can write 



xn\ = |(i-| 



1 + 1 + 0(0 



(89) 



that has the same e-expansion as Eq. (83), but it is expected to converge more rapidly to the 
correct result. From Eq. (89) we get for the two-dimensional Ising model -X"|?[l] ~ 0.20, that 
is much lower than the value that has been determined so far X^ ~ 0.33 [26,29,11,30,32]. 
However, we stress that a robust field-theoretical prediction for X^ for the two-dimensional 
Ising model requires a (difficult) higher-loop computation. For the three-dimensional Ising 
model we obtain -X"|?[l] ~ 0.37, to be compared with the direct estimate -X"|?[u] ~ 0.39. 
Note that, as usual in d = 3, the spreading of the different estimates is much smaller, 
signaling a higher reliability of these predictions. (We recall that, in three dimensions, the 
field-theoretical estimate for the FDR of the magnetization is X^[21oop] = 0.429(6) [26].) 

In Fig. 5 we report Eq. (89) for < e < 3 compared with the direct estimate Eq. (83) 
and with the two- loop FDR of the magnetization [26] . 

C. The 0(N) model with N > 2. 

Eq. (83) allows us to provide predictions for the O(N) models with arbitrary N. So 
far the non-equilibrium dynamics of models with continuous symmetry (N > 1) has been 
numerically studied only for the XY (N = 2) model, both in d — 2 [21,31] and d — 3 
[33]. The value that has been determined in d — 3, i.e., X^ = 0.43(4) [33] is in good 
agreement with the two-loop field-theoretical prediction X^ = 0.416(8) [26]. The observable 
dependence of Xq has not yet been addressed in these cases. 

To obtain an estimate of X^ T for the three-dimensional O(N) model we again constrain 
the e-expansion at the lower-critical dimension (d\, c ,d. = 2 in this case). We assume X^(d = 
2) = for N > 2. This surely holds for N > 2, since for d — 2 these systems are in the 
coarsening regime. Such an assumption is instead questionable for N = 2, where it is known 
that a finite-temperature phase transition of topological nature takes place at Tkt- A linear 
constraint would lead to 



*yiH(i-5 



16 - N -e + OU*) 



6(^ + 8) 
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(90) 



On the other hand the exact results for the spherical model [13] suggest that the approach to 
d — 2 is quadratic rather than linear, so it is tempting to implement a quadratic constraint 
also for finite N 



A?w = i(i-f s 



20 + N -e + o(S 



3(iV + 8) 



(91) 



In Fig. 5, we report Eqs. (90) and (91) for < e < 2 and N = 2, 3, oo. It is apparent that 
for N = oo the approximation obtained implementing the quadratic constraint reproduces 
better the result for the spherical model. For the three-dimensional XY model (N = 2) we 
get X|? ~ 0.37 from direct estimate (no constraint, [u]), X|? ~ 0.31 from linear constraint 
[1], and X^? ~ 0.22 from the quadratic one [q]. For the three-dimensional Heisenberg model 
(JV = 3) we find, instead, Xf [u] ~ 0.35, X%[1] ~ 0.30, and Xf[q\ ~ 0.21. 

Even the results with constraints are rather scattered, making very difficult to provide 
robust estimates in d = 3. However, we can surely conclude that X^ < X^ for all 2 < d < 4 
and that their difference should be large enough to be observed in Monte Carlo simulations 
in three dimensions. The analysis of the non-equilibrium behavior within the e = d — 2 
expansion [45] may clarify which, between the linear and the quadratic constraint, is the 
proper one close to d — 2, even if the results for the spherical model strongly suggest the 
latter. 

For Xj? we note that it is very close to Xjfi, even for e C 1 (see Fig. 5), making the 
numerical detection of such a difference probably very difficult. 

VI. CONCLUSIONS 

In this paper we considered the problem of the definition of a unique effective tem- 
perature via the long-time limit of the fluctuation-dissipation relation for critical systems 
quenched from a high-temperature phase to the critical point and evolving according to a 
purely dissipative dynamics. Within the field-theoretical approach to non-equilibrium crit- 
ical dynamics and by means of appropriate RG equations, we obtained the general scaling 
forms for the response and correlation functions of a generic local observable 0(t), from 
which it is possible to derive the FDR X (t, s). 

We found that in the Gaussian approximation all the local operators have X@ = 1/2, 
allowing for a definition of a unique effective temperature. This equality is broken already at 
the first order in the e expansion for the quadratic operators we considered (namely the total 
energy and the tensor, see Eqs. (41) and (42)). Let us point out that our results go further 
than those obtained for the one-dimensional Ising model [11,12] and for the spherical model 
[13]. In these cases Xq = X^ for all the observables O, except for the total energy. This 
operator is conjugated to the temperature of the bath but not to the actual one (if any) of 
the system. On this basis one could doubt that the energy operator is not as suited as others 
to define the effective temperature, resulting in a different X°°. Nevertheless we find that 
there is at least one further operator, namely T^-, having Xq? ^ X^,Xjf. This explicitly 
shows that, at variance with what is often conjectured, a unique effective temperature can 
not be defined for this class of models. 
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Our results for N — > oo are always in agreement with the recent ones for the spherical 
model [13], calling for a proof of a possible correspondence (or for a counterexample). They 
instead disagree (at qualitative level) with the available numerical simulation of the Ising 
model [11] apparently giving X^? = X^- Probably a more accurate measure of X|? is 
required to detect the difference (if any) between X|? and X^. We also provided theoretical 
predictions for the O(N) model for arbitrary N in 2 < d < 4. It should be possible to check 
them quantitatively in d — 3, where the e-expansion is expected to be more accurate. This 
calls for numerical simulations or real experiments in three-dimensional systems. 

Let us finally comment that it would be interesting to understand how the standard 
scenario of effective temperature [1,62,7] can be generalized to the case when each sector of 
a theory has a different T c g, as our results indicate to be the case for critical systems beyond 
the mean-field approximation. 
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APPENDIX A: TWO-LOOP FEYNMAN DIAGRAMS 

In this appendix we report all the details of the evaluation of two-loop integrals. Here 
and in the following we will denote £< = min{£i,£ 2 } and t> = max{£i,£ 2 }. For the diagram 
G\ one finds 



rti 

/o 
rti 

/0 

+ f^dt'lhih - - 2/ 2 (£ 2 )][/ 1 (£ 1 - t') - hit,)} 
Jo 

rti 



Gi(£i,£ 2 ) = rdt'B^t^Bcdt'M) 
Jo 

= /Vj 2 (t / + t 2 )[/ 1 (t 1 -t')-/l(tl)] 

Jo 



+e(h - 1 2 ) f W[/ 2 (t' - 1 2 ) - 2/ 2 (t')][/i(ti - - Ji(ti)] . (ai) 

Jt 2 

To isolate the dimensional poles of this expression one has to keep in mind that Ii(t) diverges 
for t — > and d — 4. These singularities have to be removed by proper subtractions. The 
remaining part can be expanded in a regular power series of e = 4 — d. Then one finds that 



Gi(t 1 ,t 2 )=/ 2 (l)/ 1 (l)f 1 |^W 



•+r-(-r 
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[(1 - x)- l+e ' 2 - 1] 



/2 \ / £ 2 \ e/2 

+e(h -t 2 ) dx -± 

Jo \t\ 



X 



e/2 ~(t 2 yi 



[(1 - x)~^ 2 - 1] 



+0(ti - t 2 ) / dx 
Jtl/tl 



t 2 Y' 2 
x 

h 



- 2x^-(l-^ e/2 



h 



[(1 - x)- 1+e ' 2 - 1] 



+0(*i - t 2 

+6(t 2 - ti) / dx 
Jo 



1- — 



UV' 2 



-2 



2 A _M £/2 
e V ti, 



t 2 
1-- 

h 



X 



£ / 2 /£ 



-0(*2 - il) 



e/2 



e/2 , tV /Z 



[(1 - x)~ 1+e ' 2 - 1] 



Expanding this expression in power of e, one gets 



(A2) 



Nil 



ti 



N 2 



ti 



Gi(«i,*2) = "f-log f 1 - f J - -^(7* + log2 - - + logt>) log f 1 - f 



+ 



^3*2 



i 



log 1 



log 1 



> 



> 



Nl 



t 



t 



t 



1 



t 



< 



l0 ^'° g [ l -i + 2 l0g 1 + i; + 2 l0g T-e 



~e(t 1 -t 2 ) 



^-T^ 2 { l -d- lo 4 lo < 1+t i 



to 



U\ . to 
— lo 

t 2 /*l 



+2Li2ll __ j+L12 1- 



Li, 



* 2 /«i 



1 + t 2 /*i 



^ 



-r^(t 2 - h 



IX 



2 Li 2 

6 



ti/h 



LU 



0(e). (A3) 



i-h/hj yi + hfa, 

For the following computations it is useful to introduce the function W(x) defined as 



rr .. . , 2 / 1\ , 1 + x 



J¥ (x) = -2 + — - 2x + 2 log(l - x 2 ) + x + - lo , 

2 V x J 1 — x 

— 21ogxlog(l — x) + log X log (1 + x) — 2 log 2 (l — x) 



-3 log 2 (l + x) - 2Li 2 (l - x) - Li, 



1 



1 + x 



Li 2 



x 



1 + xJ ' 



(A4) 



whose expansion for small x is W(x) = — 17x 2 /6 + 0(x 3 ). W(x) enters the expression of 
Gi(t, s) + Gi(s,t) as 



G 1 (t,s) + G 1 (s,t) 

N 2 1 



N 2 



e log M - — 1 _^( 7£ + log2 + logt)logll-^ 



Nl 



W(s/t) . (A5) 
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The expression of the diagram G 2 is 

G 2 (t,s)= fdt'B R (t } t')B K (t\s), 

•Is 

after the proper subtractions one finds 



(A6) 



G 2 {t,s) = li{l)t 



- y+ ' ' / dx 

ls/t 



'l-x)- l+ i - [1 



S \ - 1+ 2 



s \ - 1+ f 



x — 



-!+f 



-l+£ 



4 / sV /2 2 



'^(l-x)- 1 ^/ 2 ^- 14 - 6 / 2 -!) 



s/t 

sX e/2' 

1 



(A7) 



whose expansion is given by 



+ 0(6) . (A8) 



In order to compute the diagrams Lj one has to determine the tadpole (loop of the correlation 
function), given by 



P(t) 



d d q 



C°Jt,t) = -h(t) 



(A9) 



where the dimensional regularization has been used. The contribution of the tadpole to the 
two-point correlation function of the order parameter (i.e., (<£>(q, ti)(p(— q, t 2 )), a subdiagram 
of Li), is given by 



Dnc(ti,t 2 ;q) = J^ 1 dfffP tl (t 1 ,lf)P(t , )C^lf,t 2 ) 



e -q 2 (t 1 +t 2 -2t') 



"dt'hif 



-e(ti-t 2 )/i(l)Cj(ti,t 2 



t 2-d/2 e _ g 2 (tl+t2) ' 

/l(1) 2-rf/2 g2 



,2-d/2 ,2-d/2 
H — f-2 



2 - d/2 
While the contribution to the response function (y?(q, £)</?(— q, s)) is given by 

,t +2-d/2 _ 2-d/2 

D*R(t,s',q) = Jdt'Rlit^P^Rlit^s) = -I 1 (l)—^———R q (t,s) 

Using the previous expressions, the diagram Li can be written as 

r d d a +2-d/2 _ s 2~d/2 
L i(*, s ) = IttAzC^ s)D RR (t, s; q) = -h(l) ^— ^ B R (t, s) 



(A10) 



(All) 



(27T) 



2 -d/2 



4.2-d/2 _ Q 2-d/2 

-A 2 (i) 2 _/ /2 [{t-sf-^-e-^] 



(A12) 
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whose expansion is 



LJt, s) = ^-- — log - + 0(e) . 
v ' ; 16 tt-s 6 t w 



(A13) 



Analogously 



2-d/2 



Ht, s) = jjAaK^ s ) D Rc(s, t-q) = - f Q dt'h{t')h{t - if) + 1,(1)1^)^—— , (A14) 
and subtracting from the integrand its singular behavior for t' — > 0, we find 



L 2 (t, s) = - t'dlfh(i!)[hU - if) - hU)} = -I?(l)t 3 - d [ S,t dxx l - d/2 [{l - xf- d/2 - 1] 

JO JO 



whose expansion is 



£ 2 (M) = ^log['i s 



+ 0(6 



(A15) 



(A16) 



The expression for L3 is given by 
L 3 (t, s) = J-nA~ d K^ s ) D Rc(t, s; q) 



2-d/2 f2-d/2 _ „2-d/2 

o dt'h{t')h{t - if) + JiWiW^Td/S ~ h{l) 2-6.12 Mt > S) ' (A17) 



that, keeping into account Eqs. (A12) and (A14), leads to 

L 3 {t,s) = L 2 {t,s) + L 1 (t,s). 
The last diagram is L 4 , for which we find 

dd< i ^0 



(A18) 



Hh, h) = jj^y C l^ t2)D RC (t 1 ,t 2 ; q) 



^dt'h{t')[I 2 {t > -t')-I 2 {ti + t2-t')} 



-h(l) 



,2-d/2 



2-d/2 



[/ 2 (t>) - h{h + h)} 



,2-d/2 ,2-d/2 

-0(t 1 -t2)h(l) 1 2 ~ d J 2 B c ( tl ,t 2 ). 
Subtracting the singular part of the integrands one finds 

Hh,t 2 ) = -/ 1 (i)/ 2 (i) 



(A19) 



x 



J £ e > / c/rraT 1 ^ 2 



(1-rr)^ 2 - M + *<-x 



e/2 



- 1 






6/2' 



+e(h - 1 2 )- (tf - tf) (ti - t 2 ) e / 2 + (*! + 1 2 ) £ / 2 - 2ti 



e/2 



(A20) 
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whose expansion is 



L4[ti, ti 



m 



Li 2 



t</h 



i + tju 



Li 2 



9(h - h 



N, 






''loR-lo.d] -p 



Of 



(A21) 



Inserting the previous expressions into Eqs. (74) and (75) we end up for the response 
function with 



*•<*.»>= ii- t§K^ 



X 



e C 2a #\ i /, 



1 j_ C 2b9, t 

1 + ^4 l ° g S. 



F R (s/t) + 0(~g 2 ,e~g,e 2 ), (A22) 



and for the correlation function 

C 2a # 



C B (M) 



1 



Ci e 



^A 



s 2 1 



C' 



t t- s 



X 



i h- 1 1 - f § ) M* - .; 



1 + ~4 l ° g S. 



F c (s/t)+0(~g 2 ,e~g,e 2 ), (A23) 



where Ar, Fr, Ac, and Fc are the expressions reported in the text [Eqs. (76), (77), (78), 
and (79)]. To shorten the formulae we introduce the following functions: 



Wi{x) = log 2 (l + x) - log 2 (l -x) -2 log(l - x) log(l + x) + 2 log(l - x 2 



(A24) 



W 2 (a;) = W(x) + log 2 (l + x) + 3 log 2 (l - x) + 2 log(l - x) log(l + x) - - log(l - x 2 ) , (A25) 
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Wafa; 



(x) = Li 2 f — — j - Li 2 (x) - log(l - x 2 ) . 



(A26) 



The bare expressions (A22) and (A23) have to be renormalized according to [using the 
fact that O b = Z R , 6 B = Z d O R , tt B = (ZJZ^Qr = Q R + O (~g 2 , eg , e 2 ) , and Z d = 

{Z^Z^ 2 Z = Z + OCg 2 rge,e 2 ) = (l - gff ) + 0(^ 2 ^e, e 2 ) [45]] 

C a (t, S ) = (l + ^|)C B (t,s), 

i2 fl (M) = (l + ^f)i*B(M). (A27) 

Exponentiating the logarithms, we recover the expected scaling forms and exponents with 

A-N 



a 



C2ajf_ 

d 2 



E: 



T : 



2(N + 8 

N + A 
2(N + 8 



-e + 0(e 2 ), 
-e + 0(e 2 ), 



(A28) 



and 
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C 2b g* N + 2e 



0(e 2 ) for E,T 



(A29) 



d 4 iV + 82 
in agreement with Eqs. (50), (51) and the scaling forms (39) and (40), both for E and T. 
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